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Abstract 


Simulated equivalent width data of weak lines in occultation spectra have 
been analysed to retrieve both the atmospheric temperature profile and the 
tangent heights of the rays. It was assumed that the siixing-ratio of the 
absorbing gas was known but it is not necessary to know the satellite height. 
By analysing lines with temperature-dependent and temperature-independent 
intensities temperatures accurate to 1 I and tangent heights to 0.25 km were 
determined over the range 20 to 90 km. 
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Introduction 

There It considerable Interest In the ramot* sensing of ataospherlc 
temperature profiles, for example, temperature and pressure profiles are 
required to obtain vertical mlxlng-ratlo profiles of absorbing gases froa 
solar occultatlon absorption spectra. 

There have been only a few attempts to obtain teaperature profiles froa 

absorption spectra. Vatson and Tin (II first propored the use of C0 ? lines 

In solar spectra collected with a balloon-borne lnstruaent. Toth (2) 

developed an Iterative technique using the equivalent widths of high J lines 

of the v band of CO . Rlnsland et al. (3) have applied the spectral - 
3 2 

curve-fitting technique of Chang and Shaw (4) to lines of the 10.4 iai C0 ? 
band In solar spectra. Park et al. (5) and Park (6) have Investigated the 
simultaneous estimation of temperature and pressures from frequency-averaged 
transmittance values of C0 ? near 2385 cm 1 and 2255 cm *. 

VI th the exception of the methods described by Park (5,6) It has been 
universally assumed that all other required quantities Including the height of 
the observing Instrument and the pressure profile are known. This Is not 
necessarily true for satellite observations and It Is desirable to develop 
methods for retrieving Information about the atmospheric paths traversed from 
the experimental data. 

Previous Investigators have noted that the sensitivity of the temperature 
retrieval depends on the lower state energy (K") of the line. According to 
Gray and Selvldge (7) the Intensity S(T) of lines of linear molecules can be 
written 

S(T) - S(T )(T /T) 2 exp[l .439 B"(T-T )/TT ], 
o o o o 


( 1 ) 
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where S(T q ) Is the Intensity at a reference teaperature T q . The high J 
lines of the v 3 C0 2 R branch, which are often proposed for tesq>erature 
retrieval (2,3,5) have I" values of 1500 cm -1 or more. The fractional 
change in intensity AS/S per Kelvin for these lines increases monotonically 
free 3% at 300 K to 8% at 180 K (5). Lines with saaller lower-state energies 
are less testperature sensitive and have been used for pressure retrieval (5). 
Shaffer and Shaw (8) have suggested that these lines can be used to determine 
absorber aeounts along a particular ray. The instrument llne-of-slght can 
then be determined without requiring independent information of doubtful 
accuracy (8) provided the mixing-ratio of the absorber is known . In 
particular, the position of the satellite along the ray path is not required 
although calculations are usually carried out by assuming some arbitrary 
height and retrieving the corresponding set of zenith angles of the arriving 
rayr. for that position. Similarly it is assumed the latitude and hence thn 
radius of the earth are known and the tangent height is obtained by assuming a 
surface pressure. The implications of these assumptions have been discussed 
by Shaffer et al. (8) who have shown that for earth geometries, they do not 
detract significantly from the accuracy. 

Temperature Profile Retrieval 

We have retrieved a vertical temperature profile from the simulated 
equivale ' widths of lines with large E" values and simultaneously obtained 
the apparent zenith angles of the rays which produced the spectra by including 
the equivalent widths of additional lines with small B". 

The analysis is similar to that used by Shaffer and Shaw (9) to determine 
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the resolution of alxing-ratlo profiles. Provided a line is weak the 
contribution to the total equivalent width V froa each layer traversed by a 
ray with apparent senith angle s is additive and 


V(b) 



Sd^) dtTj.iy a^d). 


( 2 ) 


where the nixlng-ratlo of the absorbing gas p Q is assuaed to be constant 

with height, d(T ,A ) is the atsnspherlc density of the Jth layer with 
J J 

teaperature T and altitude A and At (a) is the path eleaent through the 

J J J 

layer. 

A total of eighty equivalent widths for the four pairs of lines in Table 1 
was calculated by using Rqs. (1) and (2) and dividing the U.S. Standard 
Ataosphere (10) froa 1 to 100 ka into 200 layers. For purposes of calculation 
the observer was at 500 ka and the 40 ray paths had tangent heights spaced 
about 1.7 ka apart between 23 and 83 ka. It was assuaed that useful equivalent 
widths for each pair of lines could be obtained over about 2 scale heights or 
16 ka. As indicated in Fig. 1 the ranges of tangent heights of tn« rays 
producing each pair of lines did not overlap. The line intensities were 
adjusted so that the largest equivalent width was the saae for every line. 

The quantities d(T^,A^) a^d Al^(x) in Kq. (2) were calculated with a coaputer 
prograa obtained froa Snider and Goldman (11). Gaussian "noise" was added to 
produce slaulated experimental data with a SNR of approxiaately 100 for the 
largest equivalent width values. 

Because the line intensity is a nonlinear function of the teaperature, an 
iterative technique was used to estlaate the teaperature profile and the zenith 
angles. The alnialzatlon portion of the coaputer prograa mis written by 
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Ralston and Jannrlch (12) and Is available under the nam BHDPAR (13). The 
refresslon analysis produces estimates of the parameters and their 
uncertainties by minimizing the sue of the squares of the differences between 
the experimental equivalent widths and those obtained from the parameter 
estiaates . 

In principle, eighty temperature and xenlth angle values can be obtained 
by analysing the eighty equivalent width values. However, since the 
uncertainty increases as the resolution Increases (9) only ten tesperature 
values were retrieved. The temperatures of the other layers in Kq. (1) were 
obtained by linear interpolation. 

If each spectrum is taken sequentially in time and the relative times of 
observation are known, then the senith angle corresponding to each spectrum 
can be interpolated from a small number of values, provided the zenith angle 
changes smoothly with time. These assumptions depend on the guiding accuracy 
of the optical system. In this study eleven zenith angles were retrieved and 
the remainder were calculated by interpolation . 

Initial parameter values must be supplied in iterative analytical 
techniques. The ability to obtain a solution was tested by using an isothermal 
(240 K) atmosphere as the initial guess. The initial guesses for the zenith 
angles were constructed by adding a random uncertainty to the true values 
corresponding to an average tangent height uncertainty of about 2.5 km. 

If the SNR is large and the initial parameter guesses are good a solution 
is obtained in a few Iterations. Because of the poor initial guesses in the 
case described above as many as fifty iterations were necessary before the rms 
value of the residuals corresponded to a SNR of 100 . 
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Sets of simultaneously retrieved temperatures and zenith angles are given 
in Tables 2 and 3. These temperatures are compared with those in the Standard 
Atmosphere in Fig. 2. Estimated standard deviations are given in parentheses. 
The number of temperature paraswters is the smallest which allows the principal 
features in the profile between 20 and 85 km to be reproduced. The differences 
between the "experimental** equivalent widths and those calculated from the 
retrieved parameters are uhown in Fig. 3. These differences are essentially 
independent of height, indicating that a satisfactory retrieval has been 
obtained. 

All the temperatures except those above 90 km were retrieved with an 
accuracy of about 1 K and the zenith angles to better than 0.005°, which 
corresponds to a tangent height uncertainty of about 0.25 km. It is seen from 
Fig. 1 that the equivalent width of each line decreases nearly exponentially 
with tangent height and this causes the SNR to decrease in a similar fashion. 
This accounts, in part, for the poor agreeswnt between temperatures above 
85 km. In addition, the accuracy of the estimates is reduced by the lack of 
equivalent widths corresponding to rays above 85 km. 

It was found that the correlations between the retrieved values of the 
zenith angles and the temperatures were small. This is due to the very 
different dependences of the absorber amounts along the paths with tangent 
height and the change of temperature with height in the atmosphere. However, 
the retrieved tangent heights of the rays do depend on the temperature 
profile. If the temperature profile is assumed to be known then only the 
tangent heights of the rays have to be retrieved. The effects of assuming a 
ten parameter temperature profile with incorrect values for the temperature 
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parameter at 33 ka on tha retrieved tangent heights are shown in Pig. 4. The 
greatest error occurs near 33 las. These results Indicate that accurate 
temperature profiles are required If accurate rap paths are to be determined. 

Other retrievals have been made by estimating temperatures at equally 
spaced height Intervals of 10 km. As expected, the temperature profile is not 
reproduced as well as In Pig. 2 but the differences can be attributed to the 
poor model used. 

Conclusions 

It has been shown that the temperature profile and ray paths can be 
simultaneously estimated from the slwlated equivalent widths of weak 
absorption lines of a gas whose mlxlng-ratlo Is known. The estimates depend 
on the characteristics of the spectra available for analysis including the 
number of lines, their intensities and lower state energies and the model used 
for the temperature profile. The results described here were obtained with a 
minimal number of lines. Although it has been assumed that the mixing-ratio 
of the absorber is constant with height, as for example CO^, or 0^, it is 
only necessary that the mixing-ratio be known. 
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TABLI 1 


Intensities end lower state energies of line pairs 
used and the tangent height range for each pair. 


Line 

Pair 

Line 

Intensity (296 K) 
[csrVaol or 2 ] 

g- [<*-!] 

Tangent Height 
Range [la] 

1 

7.2 z 10” 22 

2000 

22 * 38 


3.9 X 10" 23 

0 


2 

2.3 X 10" 21 

2000 

40 ■* 34 


3.3 X 10“ 22 

0 


3 

2.3 x 10" 20 

2000 

36 ♦ 70 


3.7 x 10“ 21 

0 


4 

1.2 * 10- 1# 


71 89 


2.3 * 10" 20 

■1 
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1ABLC ? 


Ntrl«v«d Temperature* (K) Compared with Corresponding Value* 
of the U.S. Standard Ataospfcere* Temperature Prof 11a 


ratur* 


Altitude 

(la) 


Initial 

Quess 


Retrieved 

*r 


Standard 
Atmosphere T c 


Difference 
(Tp - T C > 


100 

95 

86 

70 

60 

61 

48 

40 

33 

20 


240 

240 

240 

240 

240 

240 

240 

240 

240 

240 


246 ( 280 ) 

150 ( 20 ) 

186 . 8 ( 1 . 3 ) 

218 . 0 ( 4 ) 

247 . 6 ( 4 ) 

271 . 5 ( 9 ) 

271 . 8 ( 6 ) 

250 . 3 ( 4 ) 

230 . 1 ( 4 ) 

216 . 2 ( 6 ) 


196.1 

188.4 
186.9 

219.6 

247.0 

270.6 

270.6 

260.4 

231.0 

216.6 


60.6 

- 38.4 

- 0.1 

- 1.6 

♦ 0.5 

♦ 0.9 

♦ 1.2 

- 0.1 

- 0.9 

- 0.4 










TABLE 3 

Tha alavan astisatad smith anglas usad to aodal 
tha AO ray paths of tha sin a la tad data. 


Zmith Anglas (Dagraas) 

Tangant 

1 I I Haights 


Initial Ratriavad Corract Diffaranca (la) 

Ouass s R s c (sp - s c ) 


109.95 109.904(1*) 110.000 -0.016 85 
110.196 110.153(4) lit 149 -0.004 79 
110.241 110.299(2) 110.k.«O 0.000 73 
110.387 110.448(2) 110.448 0.000 67 
110.532 110.597(1) 110.597 0.000 60 
110.778 110.747(2) 110.747 0.000 Ml 
110.786 110.893(2) 110.896 -0.003 48 
111.069 111.047(1) 111.045 *0.002 42 
111.115 111.200(3) 111.195 *0.005 35 
111.360 111.342(1) 111.344 -0.002 29 
111.506 111.495(1) 111.493 *0.002 24 
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figure Caption* 

PIG. 1 The simulated equivalent width data calculated from Kq*. 1 and 2 using 
the line parameters in Table 1 and the U.S Standard Atmosphere. 

PIG 2 A. U.S. Standard Atmosphere temperature profile, showing initial 

240 )( guess of the tesqwrature profile. B. Differences between U.S. 
Standard Atmosphere Temperature profile and temperature calculated 
from parasmter estimates in Table 2. 

PIG 3 Residuals between the estimated widths V calculated using the U.S. 

© 

Standard Atmosphere Temperature profile and correct raj paths and the 

widths V calculated using the retrieved temperature and senlth 
c 

angle parameters in Table 2 and 3. 

PIG 4 Changes in the senith angle retrieved with the s»del temperature 

profile fixed, Including a systematic error AT in the parameter of 


altitude 33 km. 



Tangent Height [km] 
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Residual (cje-<j c ) [cm' 1 ] 


F.j. 7 





A Tongent Height [km] 




